Draft version June 15, 2012 

Preprint typeset using I^^'T^]X style emulateapj v. 5/2/11 



MODELING THE RED SEQUENCE: HIERARCHICAL GROWTH YET SLOW LUMINOSITY EVOLUTION 
Rosalind E. Skelton,^ Eric F. Bell^ and Rachel S. Somerville ^ 

Draft version June 15, 2012 

ABSTRACT 

We explore the effects of mergers on the evolution of massive early-type galaxies by modeling the 
evolution of their stellar populations in a hierarchical context. We investigate how a realistic red 
sequence population set up by z ^ 1 evolves under different assumptions for the merger and star 
formation histories, comparing changes in color, luminosity, and mass. The purely passive fading of 
existing red sequence galaxies, with no further mergers or star formation, results in dramatic changes 
at the bright end of the luminosity function and color-magnitude relation. Without mergers there is 
too much evolution in luminosity at a fixed space density compared to observations. The change in 
color and magnitude at a fixed mass resembles that of a passively evolving population that formed 
relatively recently, at z ^ 2. Mergers amongst the red sequence population ("dry mergers") occurring 
after z — \ build up mass, counteracting the fading of the existing stellar populations to give smaller 
changes in both color and luminosity for massive galaxies. By allowing some galaxies to migrate from 
the blue cloud onto the red sequence after z = 1 through gas-rich mergers, younger stellar populations 
are added to the red sequence. This manifestation of the progenitor bias increases the scatter in age 
and results in even smaller changes in color and luminosity between z = \ and z = at a fixed mass. 
The resultant evolution appears much slower, resembling the passive evolution of a population that 
formed at high redshift (z '-^ 3 — 5), and is in closer agreement with observations. We conclude that 
measurements of the luminosity and color evolution alone are not sufficient to distinguish between 
the purely passive evolution of an old population and cosmologically motivated hierarchical growth, 
although these scenarios have very different implications for the mass growth of early-type galaxies 
over the last half of cosmic history. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: general — 
galaxies: interactions — galaxies: luminosity function, mass function 



1. INTRODUCTION 



Galaxies in the local universe are found in a bimodal 
dis tribution in color- mass (or magnitud e ) space (Strat- 
eya et al. 2001 ; Kauffn rann et al. |2003[ [Blanton et al 



2003, ,Baldry et al.,,2004| . The so-called blue cloud is 
made up mostly of star-torming, disk galaxies, while the 
majority of quiescent galaxies with early- type morpholo- 
gies lie on a tight relation between color and mass known 
as the red sequence. The tilt in the color-mass relation 
(CMR) of early-type galaxies can largely be attributed 
to an increase of metallicity with m ass, with more metal 



rich g alaxies having redde r colors (Kodama & Arimoto 



19971 Gallazzi et al.||2006). The distinction between red 



sequence a nd blue c ioud~galaxies is particularly c lear at 
< 1 (e.g., [Bell e tal. 2004; Borch et al. 2006; Faber et al. 



2007 Brown et al. 2007p7and using optical-NlR colors 



or mid-lR diagnostics the populations can be separated 
reasonably well out to z ^ 3 (Labbc et al. 2005; Williarosl 

jiijr 

cloud 



et al. 



2009| [Brammer et al.|2009j |Whitaker et al.||201l| ; 
lOUgh most star tormat ion occurs in blue cloui 



iOt 

galaxies, the total mass in galaxies on the red sequence 
has grown by an order of magnitude since z ^ 2 ( Arnouts 



by a factor of t wo since z 
this 



et 



1 (the epoch we focus on in 

is paper- e g.,|Bell et al.|2004lp[ber et al.|2007| |B 
al.|2007[ ). [Bell et al.| ( |2007p showed that the increa 



rown| 

increase in 



mass on the red sequence can largely be explained by an 
infiux of blue cloud galaxies t hat have recently had their 
star formation shut down. Barker et al.|(2006) and |Ruh- 



land et al. ( |2009) used stellar population synthesis mod- 



els to show that the scatter in a red sequence that is built 
up through the quenching of star formation in blue cloud 
galaxies is consistent with observations of the evolution 
of the CMR. These works did not make any assumptions 
on what shuts the star formation down. The merging of 
gas-rich disk-li ke galaxies is a likely mechanism for rapid 



morphological (Toomre 1977 Schweizer & Seitzer 1992_ 
Barnes fc Hernquist|199D ) and color transformation, with 



some kind of feedback process usually invoked to prevent 
further gas cooling and star formation (e.g., act ive g alac- 
tic nucleus (AGN) "radio mode" feedback; Croton et al.| 



2006 Bower et al. 2006 



Somerville et al.j^200S|F 



^ ^ Most 

massive early-type galaxies are spheroidal m shape sug- 
gesting that they have for med through major mergers 
van der Wei et al. 2009b). Models in which spheroids 



ittllim lilbert et al.||2oTo| IBrammer et al.||2taFW form m gas-rich mergers, with a period of related quasar 
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activity, have been fairly su ccessful in matching a num 
ber of observation s 
20061 IBower et al 



et al.| 2008al|2009c) 



Benson et al. 2003 Croton et al 



"2006; ,SonTerville et al.p0(]8| [Hopkins 



We maintain, and will discuss below, that the growth 
in mass on the red sequence has been well measured only 
for galaxies near the knee of the luminosity function (LF; 



2 



Skelton et al. 



~ L*); it is still unclear whether the massive end of the 
red sequence has seen significant growth since z ^ 1. 
If such growth occurs, it is likely to take place through 
the merging of gas-poor galaxies that already lie on the 
red sequence since there are relatively few luminous blue 
galaxies fro m which the mo st massive early-type galaxies 
could form (Blanton 2006). The role of mergers in the 
build up of massive early-type galaxies, particularly over 
the epoch from z = 1 to 2: = 0, is strongly debated. 
In this paper we look at how a red sequence in place 
by z ~ 1 evolves with and without further mergers and 
star formation. We show that observations of color and 
luminosity evolution that have been taken as support for 
a limited role for mergers in the past are consistent with 

the changes we find using hierarchical models. 

The outline of the paper is as follows: Section [O] dis- 
cusses previous work on the evolution of early-type galax- 
ies, the support for a hierarchical picture of early-type 
evolution, and why early-type galaxies were thought to 
evolve passively from high redshift. Section [2] describes 
the mode l set up, building on the simple toy model pre- 
sented in Skelton et al. (2009). In Section Is] we present 
the resultant evolution of the color, luminosTcy and mass 
distributions in the models and compare to observations. 
We consider the implications for massive galaxies in par- 
ticular and compare the changes in color and magnitude 
to those of an ancient simple stellar population. In Sec- 
tion[5]we discuss how the various model choices affect our 
results and compare to previous work. The conclusions 
that can be drawn on early-type evolution are summa- 
rized in Section|6] We adopt a cosmology with — 0.3, 
flA = 0.7 and Ho = lOO/i km s'^ Mpc'^ with h = 0.7 
and use Vega magnitudes throughout. 

1.1. Background 

Gas-poor or "dry" mergers have been observed both 
in the local universe ( van Dokkum|2005 Mcintosh et al. 
"2008; Tal et al.||2009D and at higher redshifts (iBell et al" 
^6bi |Lotz et al.||i^Li06HLin et al.||2008| | Williams et al. 
' 201ip , and a non-negligible fraction 01 massive galax- 
les are found to occur in close pairs that are likely to 



merge (Bell et al. 2006a [Robaina et al. 



2010 



Williams 

et al. 201 IL although see Masjedi et al. 2006 who mea- 
sure very low merger rates for luminous red galaxies, 
LRGs). Recent measurements place the fraction of 
M« > 5 X 1O^°M0 galaxies in close pairs at 1 to 3% 
out to z — 1.2, implying that galaxies with masses of 
> IQ^^Mq have undergone 0.7 m ergers on average since 



this time (Robaina et al. 2010). This merger fraction 
seems to be sufficient to explain the gro wth of the red 
sequence at the massiv e end since z = 1 ( Robaina et al.^ 
2010 Man et aL]|2012 ). A particularly large uncertainty 



m interpreting measurements of the merger fraction as 
a growth rate is the timescale over which mergers can 
be identified. The timescale will be different for each 
measurement technique (e.g., visual or automated mor- 
phological selection, close pairs) and is dependent o n the 
properties of t he progenitor galaxies and the orbit (Lotz 
et al.|[2010a|b ). Recent work that calibrates the observ 



ability timescale for each method using high-resolution 
simulations in combination with cosmological models has 
shown a reassuring consistency between a number of dif- 
ferent meas urements of the m ajor merger rate and model 
predictions ( Lotz et aL]|2011 1. 



Less direct approaches for estimating the extent to 
which dry mergers build up massive red galaxies, such 
as measuring their clustering and size evolution, support 
a degree of merging activity since z = 1 but are not yet 
conclusive. Massive early-type galaxies have evolved sub- 
stantially in size o ver this period (a factor of two growth 
for a given mass, van der Wei et al. 2008). This can 
be accounted for by the increase in size of individual red 
galaxies through dry mergers, combined with an increase 
in the number dens ity of the red sequenc e galaxy popula- 
tion over time (van der Wei et al.|2009a |. Minor mergers 
may be particularly important tor the size evolution of 
early-type galaxies, as well as a viable means of explain- 
ing the evolution of the yeloci t y dispersion function since 



1 ( Bezanson et al. 2009i van der Wei et al 



Maab et al.| |2009) [Hopkins et al. 2009a, 20TU|' 
ing analyses usmgXRGs in the S loan Digi tal S 



2009a 
Cluster 
ky Survey 

(SDSS) and 2SLAQ sur veys (W ake et al.i2008. ) and the 
Bootes field (White et al. 2007) showed little evolution 



in the clustering amplitude with redshift. Interpreting 
these results in the light of the halo model suggests that 
LRGs do not evolve purely passively but that a signifi- 
cant fraction of satellite galaxies (30%-50%) must either 
merge with the central LRG or be d i srupted into the in- 
traclu ster light ( Conroy et al.||200"7 ). [Tojeiro fc Percival 
( 2010 1 also exaniine the clustering evolution of LRGs and 
find evidence for luminosity growth that is inconsistent 
with passive evolution. The growth appears to occur 
preferentially at the faint end of their galaxy selection, 
supporting an inf lux of galaxies into the sa mple through 



mergers (see also Tojeiro &: Percival| 2011 ) 



To determine how the massive end of the red sequence 
evolves, one would ideally want to measure the mass 
function (MF); however, this requires combining the 
measured LF with knowledge of the mass-to-light ratio 
(M/L), which also evolves. Stellar population synthe- 
sis models are usually used to estimate M/L, requiring 
assumptions on the age and star formation histories of 
the galaxies. To measure these quantities, the evolution 
must be known — this is one of the reasons that passive 
evolution (the evolution of a simple stellar population 
formed within a short time at high redshift, with a sin- 
gle metallicity) is often assumed for early-type galaxies. 

There is still large uncertainty in both the measure- 
ment of the number density evolution and its interpreta- 
tion. If one corrects the observed LF for the evolution ex- 
pected for a purely passively fading old population, there 



z ~ 1 (Cimatti et al.l2006VWake et al.|2006 


Brown et al. 


2007 


Cool et al.|2008; Bancrji et al.|20l0 


). J 


L'hese results 



sive red galaxies has not changed substantially over the 
last 8 billion years and that much of the mass is already in 
place by this timej^ If the early-type galaxy population 
is continuously supplemented by quenched blue cloud 

* Perhaps contrary to expectation, one of the largest difficulties 
in this analysis is determining the bright end of the LF at low 
redshift to high precision. Bright galaxies in the local universe 
have extended outer envelopes (possibly caused by dry merging) 
making it difficult to estimate t he total magnitude and accurately 
subtr act the background (e.g., [Bernardi et al.|[2"007| [Lauer et al.| 
|2007| . Small photometric errors translate into a large error m 
number density because of the exponential drop of the LF toward 
the bright end. 
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galaxies, the assumption of slow luminosity evolution 
made in the interpretation the evolution of the LF will 
be incorrect. The build up of stellar mass through merg- 
ers may act to compensate faster luminosity and color 
evolution, resulting in changes that look very much like 
passive evolut ion. This is related to the so-called "pro- 
genitor bias" ( van Dokkum fc Franx|200T ); galaxies that 
have recently moved onto the red sequence through mor- 
phological transformation will be missing from early-type 
galaxy samples at high redshift, so that at early times, 
any sample of early-type galaxies consists of only the 
oldest progenitors of today's early-type ga laxies. Con- 
vertin g the LF of red galaxies measured by [Brown et al. 
(|2007[) to a MF using the M /L evolution from Funda- 



mental Plan e measurements (van Dokkum & van der 



Marel 20071 rather than using ste llar population mod - 



els gives rise to stronger evolution ( Robaina et al.|[2010 ) 



Other authors have recently found a similar mcrease m 
number density (approximately a factor of four) at the 



mass ive end using the galaxy MF (Matsuoka & Kawara 
2010|. A recent analysis that uses clustering measure 



ments to select likely progenitors of early-type galaxies 
at high redshift measures a factor of ~ 5 growth in num- 
ber dens ity, supporting subs tantial growth through dry 
mergers (Padilla et al. 2011). 

The uniformity ot the stellar populations in early-type 
galaxies suggested by the tight CMR in clusters has long 
been used as an argume nt that they formed relatively 



1997 Bower et al. 



quickly at high reds hift ( Ellis et al., , 

[l998^ Andreon||2011 ) . However, receiit work on the evo- 



lution ot the red sequence has shown that it is difficult to 



a purely passive, old model 1 


Kriek et al. 


2008 


Ruhland 


et al. 


2009 


Whitaker et al. 


2010). 'i'he sma 


1 scatter 



in the colors of early-types in clusters was thought to 
impose a limit to the amount of gro wth that could oc- 



cur vi a merging (Bower et al. 1998). In Skelton et al 
(2009 SBS09 hereafter) we explored how the amount ot 



dry merging expected in a standard hierarchical model 
affects the scatter and slope of the CMR. We assumed 
that major gas-rich mergers move galaxies from the blue 
cloud onto the red sequence, placing them onto a "cre- 
ation red sequence" given by the local relation between 
color and magnitude. We calculated where red galaxies 
at low redshift lie in color-magnitude space by combining 
the fluxes of their progenitors, assuming that merging of 
galaxies on the red sequence does not involve further star 
formation. By using the local CMR as a reference, the 
model essentially assumed that early-type galaxies are 
evolving passively and at the same rate. With this sim- 
ple model we showed that dry merging can reduce the 
scatter and flatten the bright end of the CMR. The re- 
sulting relation was found to agree well with the observed 
red sequence from the SDSS. 

In this paper we incorporate more realistic evolution of 
the stellar populations into the model described above to 
test how mergers affect the color, luminosity, and mass 
evolution of the red sequence since z — \. Our aim is not 
to make a complete, accurate model of the galaxy distri- 
bution; the point is rather to show how the progenitor 
bias expected from hierarchical growth in a standard cold 
dark matter cosmology affects the evolution of massive 
galaxies. T he models, very mu c h in the spirit of ear- 
lier work by Bower et al. (19981; van Dokkum & Franx 



( |2001| ; [Blanton] ( |2006[ ); [Skelton et aL] ( [2009[ ), and oth- 
ers, use stellar population synthesis niodeling to follow 
the luminosity evolution of galaxies, but set the model 
galaxies into a hierarchical context using merger trees 
extracted from a semi-analytic model (SAM) of galaxy 
formation. We assume that the red sequence is built up 
mainly through major gas-rich mergers. With a reason- 
ably realistic red sequence population in place at z — 1, 
we explore the consequences of subsequent gas-poor and 
gas-rich merging. We compare this to the evolution of 
the same population undergoing only passive evolution 
with no mergers after z = \. We show that at fixed 
mass, dry merging and recent additions to the red se- 
quence decrease the change in color and luminosity. The 
resultant evolution is consistent with that of a popula- 
tion that formed early and evolved passively, even though 
there has been significant merging activity. 

2. MODEL FRAMEWORK 

To investigate how the colors and magnitudes of galax- 
ies evolve in a hierarchical universe, we incorporate the 
evolution of stellar populations of different metallicities 
into the model described i n SBS 09. Galaxy merger trees 
from the [Somerville et al.[ ( [2008| S08 hereafter) SAM pro- 
vide the hierarchical base tor the models. Rather than 
using the full SAM, a complex model that has success- 
fully reproduced a number of observations, we aim to 
disentangle the effects that mergers have on early-type 
galaxy evolution with a simplified model setup. The 
SAM provides cosmologically motivated merger statistics 
that have been shown to match weh with observati ons 
( [Beh et al.|2006a[ [Jogee et al.[[2009[ [Lotz et al.[[2011 ) . 

In the S08 SAM the dark niatter merger trees are de- 
termined using the extended Pr ess-Schechter formalisrn, 
as described in Somerville et al. (!2008, 2001,) ; ,SomervilIe] 
& Kolatt ( 1999p . Analytic prescriptions for gas coolmg 
(S08, Section 2.2), star formation (S08, Section 2.5), su- 
pernova feedback (S08, Section 2.7) and AGN feedback 
from winds (SOS, Section 2.10) and radio jets (S08, Sec- 
tion 2.11) determine how galaxies form and evolve within 
the dark matter halos. A modified Chandresekhar dy- 
namical friction formula that takes into account mass 



loss due to tidal stripping ( Boylan-Kolchin et al. 20081 
is used to calculate the time taken for galaxies to merge 
after the merger of their dark matter halos. Quench- 
ing occurs when the central black hole becomes large 
enough to prevent further gas from cooling through ra- 
dio mode feedback. This often occurs for the first time 
immediately following a major merger, as these mergers 
trigger a phase of rapid accretion leading to significant 
growth of the black hole. Due to the self-regulation of 
black hole accretion, there is a correlation between the 
mass of the remnant central black hole and the mass of 
the bulge component of the galaxy, providing a natural 
link between quiescent galaxies and spheroid-dominated 
morphologies. At high redshifts and in low mass halos, 
further cold gas may be accreted via cold dense filaments 
that are impervious to heating by the radio jets, so star 
formation may begin again some time after a major wet 
merger has occurred, even if the existing gas is removed 
during the merger process. These physical processes pro- 
vide the motivation for our simple model but are not 
included in any detail. 
We construct galaxy merger trees from the SAM out- 
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Skelton et al. 



put, using the halo and galaxy identity numbers to track 
galaxies as they merge, and record the stellar, cold gas 
and dark matter masses for the two progenitor galaxies 
involved in every merger. These are used as input to 
the toy model. We infer the metallicity of each merging 
galaxy from its stellar mass using a simple approximation 
to the mass-metallicity (M-Z) relation given by 



-0.082 + 0.0098 log 



A/* 



(1) 



This relation was obtained by fitting a straight line 
through t he median metallicity as a function of stellar 
mass from Gallazzi et al. (2005) and adjusting the nor- 
malization by 0.U04 so that the color-magnitude relation 
of the model approximately matches the observations at 
z=0.9. Solar metallicity is 0.02. We assume that the M- 
Z relation does not evolve and the metallicity of each 
galaxy remains the same for the duration of its star for- 
mation. Although these are clearly simplifications, they 
are defensible choices in the light of this model frame- 
work. The existence of a mass-ga s metallici ty relation 
has been established out to z ~ 3 (Erb et al. 2006), but 
it is very difficult to pin down how the relation evolves 
observationally and there is larg e uncertainty in the cali - 
bration of the observed relation ( Kewley fc Ellison|2008 ) . 
Erb et al. ( 2006 1 find a factor oi 2 ottset between the re- 
lation at z ~ 2 and z ~ but note that the uncertainty 
is at about the same level. The rate at which the relation 
evolves seems to be mass dependent, with very little evo- 
lution since z ^ 0.8 for massive galaxies dSavaglio et al. 
2005] IZahid et al.|[20TT] but see jMoustakas et al.||20TT ' 
who find lower metallicities by 30%-60% at z = 0.7 with 
little mass dependence). 

In the toy model, the red sequence forms from the rem- 
nants of major gas-rich mergers. We begin by identifying 
the major wet mergers occurring along every branch of 
each merger tree. Major mergers are defined to have a 
mass ratio between 1:1 and 1:4, where the total baryonic 
and dark matter mass withi n twice the Navarro -Frenk- 
White (NEW) scale radius (Na varro et al.|1997 ) is used 
to determine the ratio (see S08). i'he radius of interest 
corresponds to approximately 60 kpc for a Milky Way- 
size halo. To distinguish between wet and dry mergers, 
we assume a gas fraction threshold of 20%, where the 
gas fraction is defined as the ratio of cold gas mass to to- 
tal baryonic mass (cold gas and stars). This distinction 
is somewhat arbitrary but the region of color-magnitude 
space dominated by galaxies with lower gas fractions cor- 
responds closely to the red sequence in the full S08 SAM. 
This threshold also produced the closest match between 
the model and observed red sequences in SBS09. If ei- 
ther of the galaxies has a gas fraction of more than 20%, 
the merger is classified as wet. 

We assign each galaxy involved in a major wet merger 
a star formation history, assuming that the star forma- 
tion in both progenitor galaxies was truncated at the 
time of the merger. The masses and fiuxes of the two 
progenitors are combined and evolved from the time of 
the merger to the time the remnant galaxy is "observed" . 
If further dry or minor mergers occur before the time of 
observation, the flux (mass) of the secondary galaxy are 
added to the flux (mass) of the main progenitor. The sec- 
ondary galaxy's properties are determined by following 



its merger tree back until a major wet merger occurred, 
and so on, recursively, for each branch of the merger tree. 
If no major wet merger occurs along one branch, star 
formation is assumed to have stopped at the first merger 
where the main progenitor is recorded as gas poor. Only 
mergers where the mass ratio is between 1:1 and 1:10 are 
considered. 

In keeping with the simplicity of the model, we choose 
a constant star formation history for each galaxy from 
some formation redshift z/ up until the time at which 
star formation stops, with the addition of a burst of star 
formation at the time of the merger. We present models 
with z/ = 4 below, but find that the exact choice does not 
strongly affect the results (see Section [5|. In principle, 
any star formation history could be used; an exponen- 
tially declinin g func tion is a common choice, fo r example 
(see, e.g ., iRuhland ct al. 2009). We use the Bruzual & 
Chariot (2003) stellar population synthesis models with 
a Cha brier initial mass function (IMF) and the Padova 
1994 stellar libraries to determine how the flux changes 
with time for each population of stars that form and take 
the loss of material back into the interstellar medium 
into account. We interpolate between the simple stellar 
population evolutionary tracks that are provided for six 
different metallicities to determine the appropriate evo- 
lution given each galaxy's metallicity, calculated using 
Equation ([T|. 

The star formation rate is normalized so that the to- 
tal mass in stars remaining at the time of the merger is 
equal to the stellar mass of the galaxy given by the SAM 
at that time. The additional burst mass is given by the 
cold gas mass of the merging galaxy, assuming that all 
the gas is used up in a burst lasting 1 Gyr. The nor- 
malization of the ME was found to be too low compared 
to observations if no burst was included. The gas mass 
can be a substantial fraction of the baryonic mass, par- 
ticularly at high redshifts, and for low and intermediate 
galaxies even at z < 1. We test how different star forma- 
tion histories and the addition of bursts affect the results 
in Section [Sl 

Major wet merger remnants redden and fade after the 
merger, moving onto the red sequence some time after 
star formation has stopped. In order to compare directly 
with the observe d distribution , we us e the same color cut 
that was used by Brown et al. ( 2007 ) to determine which 
of the model galaxies are on the red sequence at the time 
of observation, as described in more detail in Section [3j 
With the assumptions described above, the model red se- 
quence approximately matches both the slope and nor- 
maliza tion of the observed color-raagnit ude relation at 
z ~ 1 ( |Brown et al |2007[ [Bell et al.|2004[ see Section |3|. 
Having set up a tairly realistic population of galaxies at 
this redshift we explore how different merging scenarios 
and star formation histories affect the evolution of the 
population between z = 0.9 and z — 0.1, where obser- 
vational results are readily available for comparison. We 
present the results for three illustrative models for the 
evolution of the red sequence over the last 8 billion years, 
pictured in Figure [T] We follow the procedure outlined 
above for all three models until z ~ 1 and then allow the 
models to diverge as described below. 

2.1. Model I 
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z = 



z= 1 



z = 4 




Model I 



Model II 



Model III 



Figure 1. Schematic illustrating the evolution of a galaxy that has undergone a major merger at high redshift and a second merger at 
^ < 1 in the three models. Arrows with a color gradient indicate that the stellar populations are evolving passively. Continuing star 
formation is indicated by the blue arrows. Star formation stops at the time of the first major merger. In Model I, no mergers take place 
after z = 1, as shown by the shaded path. In Model II, star formation is quenched at z ^ 1 and galaxies that do not undergo a merger 
before z = 1 start to fade passively at this time. In Model III, such galaxies continue to form stars until they undergo a major wet merger. 
The small red arrows show approximately when the galaxies would be classified as red sequence galaxies. 



The first situation we consider, Model I, is how the red 
sequence evolves if star formation in red sequence galax- 
ies is shut off at z 1 and there are no further mergers 
of any kind. We do not allow any further additions to 
the red sequence through the shutting down of star for- 
mation in blue galaxies. The galaxies already on the red 
sequence merely evolve passively, becoming fainter and 
redder, as shown in the first panel of Figure [T] The MF 
does not evolve because the number density of galaxies 
on the red sequence is not increased through the addition 
of wet merger remnants, nor is its composition changed 
by dry mergers between existing red sequence galaxies. 
This case is extreme in the context of the standard hi- 
erarchical picture of galaxy formation, and is certainly 
unrealistic at intermediate masses, where there must be 
an influx of blue cloud galaxies after z = 1 to account 
for the observed increase in mass on the red sequenc e 



( Bell et al 
At the brig. 



2004 Faber et al. 2007 Brown et al. 2007). 



it end, which we are most interested in here 



it shows how already existing massive early-type galaxies 
would evolve without undergoing mergers. This amounts 
to the purely passive evolution often advocated for early- 
type galaxies, albeit a population in which some galaxies 
continued forming stars until z — I. 

To implement this model we truncate each merger tree 
at 2 = 1. We follow the evolution of the major wet 
merger remnants that formed before z = 1, adding the 
mass and flux contributions from subsequent mergers 
only if they also take place before z = 1. The exist- 
ing stellar populations evolve passively, with no further 
star formation, until z = 0.1. 

2.2. Model II 

In the second case (Model II) we assume that early- 
type galaxies on the red sequence by z ~ have had 
their star formation quenched by z '--^ 1, regardless of the 



mechanism. We allow all the mergers predicted by the 
SAM to occur after z = 1, but as the star formation in 
these galaxies has already been truncated, all the mergers 
occurring thereafter are essentially dry. In this model, 
which is illustrated in the second panel of Figure [l] the 
red sequence continues to grow through mergers after 
z — 1 and more massive galaxies are built up through 
dry mergers, but the red sequence galaxies we see at low 
redshift all have old populations that formed at least 8 
billion years ago. 

In Model II, major wet mergers are no longer the cause 
of the quenching of star formation after z = 1. Galaxies 
involved in mergers that have not already had their star 
formation quenched through a major wet merger by this 
time are assumed to have had a constant star formation 
rate for 4.5 Gyr from Zf — 4 to z — 0.9 and a star 
formation rate of zero thereafter. With no new source 
of fuel, this is a plausible timescale on which the fuel 
initially available for star formation could be used up, 
but the exact time at which we halt star formation is not 
critical. This particular value was chosen in order to be 
able to directly compare the changes in the three models 
from a common starting point at z = 0.9. 

2.3. Model III 

In the third model scenario. Model III, shown in the 
third panel of Figure [TJ we allow all the mergers pre- 
dicted by the SAM to occur after z — 1 and assume 
that star formation continues until a major wet merger 
occurs, rather than truncating star formation arbitrarily 
after 4.5 Gyr. This model is closest to the fufl SAM. 
Both dry and wet mergers taking place between z = 1 
and z = affect the evolution of the number density of 
red sequence galaxies. Some of the galaxies on the red 
sequence at z = will have continued to form stars un- 
til relatively recently (the time of their last major wet 
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merger) . In this way blue progenitors with young steUar 
populations are added to the red sequence throughout 
the last 8 billion years. This model illustrates the effects 
of progenitor bias. 

3. RESULTANT COLOR AND MAGNITUDE EVOLUTION 

Figure [2] shows the evolution of the CMR (top panels), 
LF (middle panels), and MF (lower panels), with each 
column representing one of the model variations. In all 
the panels of Figure [2] the dashed blue lines show the 
model galaxy distributions a,t z — 0.9 and the red solid 
lines at z = 0.1. Recall that the models are identical until 
z = I but differ in their evolution thereafter. For com- 
parison, observational results are shown as gray symbols 
and thick black lines, as described below. 

We compare the model color-magnitude relations at 
z = 0.9 and z = 0.1 with t he observed red se quence of 
galaxies in the Bootes field ( Brown et al.||2007 ) parame- 
terized by 



for a scatte r of 0.1 dex in the M-Z, as measured by I Gal 



U -V = IA- 0.08(Afy - 5 logio(/i) -I- 20) 
- 0.42(z - 0.05) + 0.07(z - 0.05)^ 



(2) 



This i s very similar to the CMR found by |Bell et aL] 
(2004) for galaxies in the Combo-17 survey. To make a 
direct comparison, we use the same magnitude and red- 
shift dependent color cut 0.25 mag below Equation d2| 
to select galaxies on the red sequence in the model. This 
is a fairly generous cut that excludes only a small frac- 
tion of model galaxies in the range of My — [—19, —24] 
in the high-redshift bin, with higher numbers of galax- 
ies excluded toward the faint end. The distributions of 
model galaxies in color-magnitude space at z = 0.9 (open 
contours) and 2; = 0.1 (filled contours) are shown in the 
top panels of Figure [2| The shaded contours represent 
1%, 5%, 10%, 25%, 50% and 75% density levels. The 
observed relations for galaxies in the Bootes field, given 
by Equation (2]), are shown as thick black lines. Al- 
though both the observed and model distributions have 
some curvature, we compare them by fitting a straight 
line to the median of the model galaxy distribution in 
magnitude bins of 0.5 mag along the relation at each red- 
shift. The scatter is calculated with the IDL routine ro- 
bust siqma using Tukey's biweight estimator (Beer s et al. 
1990 1 , which is robust to outliers (see, e.g . , Ruhla nd et al. 

Bow er et aL][T992l [Mcintosh eFatlllOOSl |Whitaker 
et al. 20101 for measurements of the scatter in the ob- 
servedTMR). 

We find good agreement between the model and the ob- 
served red sequences at z — 0.9, but note that both the 
slope and the normalization are sensitive to the choice 
of M-Z relation. A shallower M-Z relation results in 
a shallower model CMR at all redshifts. The scatter in 
the model red sequence at z = 0.9 is 0.10 mag, which is 
somewhat s maller than the observed scatter (0.21 mag) 
measured by |Ruhlandet al.| ( |2009] ) at z ~ 1 in the GEMS 
survey. The scatter decreases toward the bright end — 
for massive galaxies (M* > IO^^Mq) it is 0.08 mag. 
This is also smaller than the measurement fro m the New- 
firm Mediu m Band Survey (NMBS) at z 1 ( [Whitaker 
et al.|2010 ) . The dispersion in the CMR is largely driven 



lazzi et al.| ( |2006[ ) , increases the scatter at all redshifts, 
but still does not cause it to exceed the observed scatter. 

In the middle row of Figure [2] we show the B-band 
LFs. Observed LFs of early-type galaxies selected by 
color from t he SPSS ([Bell et al.|2003^ at z ~ 0.1 and the 
Bootes field ( Brown et al.|2007 ) at O.o < z < 1 are shown 
with gray circles and crosses, respectively. The color cuts 
we use to define the r ed seq uenc e are close to the cuts 
used by iBrown et all p007| and iBell et all (|2004|), but 



membership ot the miodel red sequence is largely deter- 
mined by merger history. As an indication of where the 
low-redshift LF lies, we show the LF of early-type galax- 
ies selecte d bv c olor from the SDSS early data release 
(iBell et al.||2 003|), corrected to Hq = 70 km s" ^ 
I'he a- band ( 



Mpc" 



he g-band ot the SDSS has been assumed to be directly 
comparable to the i?-band in the Vega system used in 
this work, and is in good agreement with t he LF estimate 
from the 2dFGRS (Madgwick et al.|2002|). Note tha t the 
observed LFs from the COM BO-17 ([Belfet^^ and 
DEEP 2 ( Faber et al. 2007 ) su rveys have lower normal- 
izations than the 'Brown et al. (2007) LF at z ^ 0.9, and 
agree better with the model at the knee of the LF. The 
bright ends are very similar in all three works. 

The model and observed MFs at z = 0.9 and z = 0.1 
are shown in the lower panels of Figure [2] The MF of 
z ~ 0.1 SDSS galaxies selected to be early- type based 
on a color cut is shown by the filled circles ( [Bell et aL 
'2003|). At high redshift (0.8 < z < 1.4) we show a recent 



by the spread in the ages of the stellar populations, as 
the model red sequence has been slowly built up by the 
quenching of star formation at different times. Allowing 



determination of the M F of red galaxies from the NMBS 
(Brammer et al. 2011[ ). This agrees very well with the 
Schechter function fit t o the MF of red seq uence galaxies 
at 0.8 < z < 1 from |Ilbert et al.| ( |2010| ). The uncer- 
tainties in the observational results for the most massive 
galaxies (Af* ^ 3 x 10^^ Mq) are large, making it difficult 
to draw a robust conclusion on their evolution. 

The model and observed number densities at the bright 
end agree fairly well at both redshifts, however the dis- 
tribution of model galaxies still differs in shape from the 
observed distribution, with too few intermediate lumi- 
nosity galaxies at low redshift. We are largely focused 
on the differences in the evolution of the bright end of 
the LF between the three models, to demonstrate the ef- 
fects of mergers, rather than trying to match the entire 
distribution of galaxies, yet it is important that the space 
density of galaxies as a function of luminosity and mass 
in the models is approximately correct at our starting 
point (z = 0.9). It was necessary to include a burst of 
star formation at the time of the merger to convert the 
remaining cold gas mass into stars, in order to match the 
observed MF. This makes little difference to the evolu- 
tion of the massive end of the red sequence after z = 1. 
We discuss the impact of bursts of star formation and 
implications of the diffe renc es between the model and 
observed LFs in Sections 15.31 and [Ol 

In Model I, which has only the passive aging of exist- 
ing stars and no mergers (left panels of Figure^, there 
is rapid evolution of the color and luminosity of galax- 
ies, resulting in a much fainter and redder population by 
z = 0. The resulting model CMR (upper left panel) is 
significantly redder than both the observed relation at 
z = 0.1 and the other models. The scatter in the rela- 
tion measured for My < —19 is small (0.05 mag). The 




Figure 2. Evolution of the red sequence from z = 0.9 (dashed blue lines) to z = 0.1 (solid red lines) in the three models. The color- 
magnitude relation is shown in the top panels, the LF in the middle panels and the MF in the lower panels. The distributions of model 
galaxies in color-magnitude space are shown by the empty contours {z = 0.9) and filled contours (z = 0.1). Linear fits to these distributions 
are compared to the observed relations fro m Brown et al. (2007), shown as thick solid lines. At z = 0.1, the LF and MF are compared to 
the observe d LFs from the SDSS (Be ll et al.|2003[ l, shown as filled circle s. At z = 0.9, the observed LF is from the Bootes field (Brown] 
|et al.|2007[ l and the MF from the INewUrm Medium Band Survey (Bram mer et al.|2011[ both shown by crosses). In the left most paneTsf 
galaxies oil the red sequence stop forming new stars at z ^ 0.9 (4.5 Gyr alter tormation) and evolve passively, with no further mergers 
occurring (Model I). In the middle panels, star formation is similarly quenched by z ~ 0.9 but mergers continue to the present day (Model 
II). Model III, where star formation may continue after z = 0.9 until a major wet merger occurs, is shown in the right most panels. 



number density of galaxies as a function of mass does not 
change after z = 1 (lower left panel). Despite the con- 
stancy of the MF, the LF undergoes significant evolution 
(middle left panel) due to the rapid fading of the stellar 
populations, which may have continued to form stars un- 
til as late as z = 0.9 in some galaxies. There are too few 
bright red galaxies remaining at low redshift, with no in- 
crease in mass through merging to counteract the fading 
of the stellar populations. This gives an upper limit for 
the amount of luminosity and color evolution; there may 
be galaxies made up of even older stellar populations 
that evolve more slowly and galaxies that deplete their 
gas and move onto the red sequence without undergoing 



mergers, which do not contribute to the red sequence in 
these models. 

In the middle column of Figure [2] the evolution of the 
model with mergers but no star formation after z = 0.9 
(Model II) is shown. Although there are still no new stars 
formed, the fading and reddening of the CMR is less dra- 
matic than in Model I because of the effects of mergers 
on both the mass and color evolution. While dry merg- 
ers increase the mass of galaxies, the slope of the CMR 
implies that remnants will be no redder than the most 
massive progenitor, reducing the slope at the bright end 
of the CMR (see SBS09). The CMR at z = 0.1 is in 
better agreement with the observed relation than Model 
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I, although there is stiU shghtly more evolution in color 
at the bright end than observed. The scatter for galaxies 
with Mv < —19 at z = 0.1 is 0.06 mag. In contrast 
to Model I, the MF of Model II evolves between z = 1 
and z — due to the redistribution of mass amongst red 
sequence galaxies and the addition of new major merger 
remnants to the red sequence. The effect of mergers is 
to shift the MF to higher masses (lower middle panel). 
The noticeable growth in the number density of massive 
galaxies caused by mergers in the models is consistent 
with the observations. The evolution of the LF incorpo- 
rates both the passive fading of the stellar populations 
and the increase in number density of massive galaxies 
from merging. These act in opposite senses so that the 
resultant change in the LF (central panel) is much less 
marked than in Model I. 

In Model III, major merger remnants added to the 
red sequence after z = 0.9 are assumed to form stars 
until the time of the merger if their progenitors have a 
high gas fraction. Such galaxies have younger luminosity- 
weighted ages and bluer colors than galaxies with stellar 
populations that all formed before z — 0.9. The result- 
ing CMR at z = 0.1 is in excellent agreement with the 
observed relation (upper right panel). The scatter in 
U — V color for galaxies in Model III with My < — 19 
at z = 0.1 is 0.09 mag, which is slightly smaller th an the 
scatter measured for SDSS galaxies at low redshift ( Ruh- 
|land et al.]]2009|). The increase in scatter compared to 
ModeTl is in line with the expectation for a steady influx 
of blue cloud galaxies onto the red sequence afte r z = 1 



dvan Dokkum fc Franx|2001[|Ruhland et al.|2009{|Harker 



et al.|2006[ ). For galaxies with > W^Mq the scatter 
in U — V color is 0.07 mag, consistent with the scatter 
measured fo r galaxies in this m ass range in clusters out 
to z = 0.84 ( |van Dokku^|2008| . 

The MF in Model 111 evolves in the same way as for 
Model II. Although some of the galaxies have a different 
star formation history, the same total mass in stars is 
added to the red sequence. Galaxies that formed stars 
after z = 0.9 may undergo further (dry) mergers once 
they are on the red sequence. In this way younger stellar 
populations are incorporated into more massive galax- 
ies. The fading and reddening of the stellar populations 
is compensated by the addition of mass from mergers, 
as well as recent additions to the red sequence. The 
changes at the bright end of the LF are smaller than ei- 
ther of the other models. The resulting evolution of the 
luminosity function is very mild, although there is sub- 
stantial growth through merging after z = 1. This gives 
an alternative explanation for the mild evolution previ- 
ously interpreted as the result of purely passive fading of 
galaxies that formed their stars at high redshifts. 

4. IMPLICATIONS FOR THE EVOLUTION OF MASSIVE 
GALAXIES 

Figure [3] summarizes the changes in color and magni- 
tude between z — 0.9 and z = 0.1 for model galaxies with 
log(M/M0) > 11. We compare the predicted evolution 
of the three models with the average changes in color and 
magnitude for passively evolving simple stellar popula- 
tions (SSPs) in the same mass range, with four choices of 
formation redshift, Zf — {2, 3, 4, 5}. To calculate the ex- 
pected evolution of the passively evolving population, we 
use the BC03 stellar population models to simulate 1000 
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Figure 3. Change in color and magnitude from z = 0.9 to z = 0.1 
for galaxies with masses > 10^ ^Mq in the models compared to the 
expected changes for passively evolving SSPs formed at redshifts 
from 2 to 5 (arrows). Model I, II, and III are shown by the open 
square, filled circle and filled triangle, respectively. 

model galaxies drawn from the observed MF for each 
formation redshift. Rather than simply assuming solar 
metallicity, we allow for the variation of metallicity with 
mass, taking into account the relative number of galaxies 
as a function of mass by assuming they are distributed 
according to a Schechter function (this is important only 
at the massive end, where the change in number density 
with mass is very rapid). We use the obs erved MF of red 
galaxies from the SDSS ( |Bell et al.||2003| ), with param- 
eters 6 



0.0107//i3Mpc-3logioM-i, logio (t^) 

10.50 and a — —0.70 to assign each galaxy a weight. This 
has been converted t o a Chabrier IMF by subtracting 0.1 
dex from the mass (Borch et al. 2006). The metallicity 
of each galaxy is calculated using Equation ([T]) , where a 
mass of log(M/M0) = 10.4 corresponds approximately 
to solar metallicity. The amount of evolution found for 
passively evolving galaxies does not depend strongly on 
the mass range (i.e., metallicity). The color evolution is 
approximately that of a solar metallicity SSP, but there 
is less evolution in the i?-band, as expected for higher 
metallicity SSPs. 

In Model I, the population of galaxies in this fixed mass 
range (log(M/MQ) > 11) is the same at both redshifts, 
whereas in the other models galaxies can move into the 
sample as their masses are increased through merging. 
The mass selection thus results in a much smaller sam- 
ple at z = 0.1 in Model I than in the other models. The 
number density of galaxies with log(M/M0) > 11 grows 
by approximately a factor of two in the models where 
merging has continued after z = 1. The population of 
galaxies is bluer and brighter at z = 0.1 than the cor- 
responding population of galaxies in Model I. Additions 
to the red sequence after z = 1 enhance the effect fur- 
ther (Model III) . The relatively small resulting changes in 
color and magnitude resemble those of a passively evolv- 
ing population that formed at higher redshift, as can be 
seen from the similarity in the position of the Model III 
point and the passive evolution vectors with z/ between 
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3 and 5. The change in color and magnitude for Model I 
replicates that of a more recently formed, passively fad- 
ing population, because there is no merging or additional 
star formation that can compensate for the rapid lumi- 
nosity and color evolution. 

Observational measurements of the change in M/Lb 
can be compared to the change in magnitude at a fixed 
mass in the models, shown in Figure[3] For cZwsier galax- 
ies with log(M^/M0) > 11 AM I Lb from z = 0.83 to 
l ies in the range of 0.8 - 1 m ag (van Dokkum et al. 



1998a| |van Dokkum fc van der Ma rcl 2007; Hold cn et a 
m)5: Wuyts ct al. 2004; va n Dokkum fc Stanford||2003[ 
van der W cl ct al. 2005 . ) , while At/ - V ~ 0.21 mag 
(van Uokkum 2008|r Most of the evolution measure- 
ments for clusters extrapolate between a single cluster at 
high z (MS 1054-03 at z = 0.83) and the Coma cluster 
at z = 0.022. The change in M/L n is larger when th e 
sample is extended to lower masses dHolden et aL][20Tol ). 
For field galaxies, there is greater spread in the measure- 
ments of d\og(M/LB)/dz from the Fund amental Plane 
( |Treu et al.|2005 [van der Wei et al.|2005[ and references 
therein). I'br Az = 0.8, as plotted in Vb'igurejsj the change 
in magnitude measured for field galaxies lies m the range^ 
of 1 - 1.9 mag (Rusin fc Koc hanck 2005; van Dokkum fc 
Ellis|2003|rvan der Wei et al.t2004,,Cebhardt et al..2003i 



2005p . Many of these studies are limited by 



Treu et a 

small sample sizes. The evolution in luminosity and color 
for all three models fall within the observed range, but 
more work is needed to pin down the observed evolution 
for a larger sample of field galaxies. 

A fairly straight-forward way of measuring the evolu- 
tion of the LF (provided that a Schechter function or 
other analytic form fits the LF reasonably well) is by 
determining the change in ma gnitude at a fixed space 



density (e.g., [Bro wn et al.|2007 |. This has the advantage 
that it is relatively insensitive to the details of red galaxy 
selection and avoids the use of a magnitude threshold 
that is sensitive to the exact shape of the LF at the bright 
end. 

In Figure [4] we show the change in magnitude for the 
three models at a space density of 10^^ Mpc^'^mag^^, 
overlaid on the observed LFs f rom the Bootes fi eld 
(NOAO Deep Wide-Fi eld Survey; [Brown et al.|[2007| ) at 
z = 0.9 and the SDSS ( |Bell et al.||2003[ ) at z =in~The 
evolution vectors for Model 1 and 111 have been offset 
slightly for clarity. The characteristic magnitudes of the 
model Schechter functions at both redshifts were shifted 
by 0.2 mag so that the model fit at z = 0.9 (thick dashed 
line) intersects with the Bootes data at the space density 
of interest. 

Figure [4] emphasizes that a red sequence galaxy popu- 
lation with a realistic mix of stellar populations at z = 1 
will undergo very rapid evolution in magnitude at fixed 
space density if there is no growth in mass (Model I, 
blue arrow). The evolution is dramatically slowed by 
growing the red sequence population through dry merg- 
ers between already-existing early-type galaxies to larger 
masses (Model II, green arrow), and slowed further by 
the addition of wet merger remnants that have recently 
had their star formation quenched (Model III, orange ar- 
row). The luminosity evolution in models that include 
mergers is reasonably consistent with the observations at 
all space densities probing the bright end of the LF. 
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Figure 4. Change in magnitude from z = 0.9 to 2 = 0.1 at a fixed 
space density of 10~'*'^Mpc~^mag~^ in the models, compared to 
the observed evolution of the LF. Th e low-z LF data (filled cir- 
cles) are from the SDSS {Bell et al.|[2CI03| ), while the high-2 LF 
data points (crosses) are from tne Uootes field l |Brown et al.|2007[ |. 
A Schechter function fit to the model LF at z = U.y, smtted by 
0.2 mag to coincide with the Bootes data at a space density of 
,f) = lO-^-^Mpc-s mag ^, is shown by the thick dashed line. The 
three models have been slightly offset in space density for clarity. 
The luminosity evolution of Model I (blue arrow) far exceeds the 
observed evolution. Models II (green arrow) and III (orange arrow) 
agree well with observations. 

5. DISCUSSION 

The behavior reported in this paper is remarkably ro- 
bust to a number of the model choices we have made. 
Though some of the galaxy properties are taken directly 
from the SAM, and thus depend on the parameters and 
prescriptions used there, others are determined indepen- 
dently in this work. Choices made in the three models 
presented here include the M-Z relation, the mass ra- 
tio used to separate major and minor mergers, the gas 
fraction threshold below which a merger is considered 
gas-rich, the gas fraction threshold above which a galaxy 
remains in the blue cloud, the time at which star for- 
mation starts, the normalization and shape of the star 
formation history, the stellar population models and as- 
sociated stellar library, and the IMF. In Model I there 
is the additional assumption of a time after which no 
mergers take place (z = 1) and for Model II, a time after 
which any ongoing star formation stops (4.5 Gyr). 

5.1. Merger histories 

We implement our model using one particular incar- 
nation of a hierarchical merger tree, taken from the S08 
SAM, but expect the results to hold more generally in 
standard hierarchical models. The merger statis tics of 
various mod els, including S ( )8, we re compared in |Jogee| 



et al.| ( |2009[ ) and |Lotz et al.| ( |2011[ ). There is broad quaT 
itative agreement between the models, but in detail the 
merger properties can differ substantially. This is most 
likely due to the different implementations of baryonic 
physics, rather than differences in the merger statistics 
of the dark matter halos, as no significant differences in 
the results of the SAM were found when the model was 
implemented within iV-body merger trees instead. 
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In our simple model, the gas fractions and mass ratios 
of merging galaxies determine which galaxies move onto 
the red sequence. The number density of red galaxies 
then depends critically on the choice of model and the pa- 
rameters used to select which galaxies are quenched, but 
we expect the relative differences between the three sce- 
narios to hold using other merger trees. It is reassuring to 
note that the baryonic mass rati o distributions of merg- 
ing galaxies in the SOS, Croton et al. (20 0^ and Stewart 
et al.l (|2009t) models are similar (Lotz et al.||201ll). The 



mean gas traction of merging gal axies m the S(J 8 SAM 



is fairly high but comparable to Stewart et al. (2009) 
The importance of mergers of diflerent gas tractions diif- 
fers from model to model as a function of time. In the 
SOS model major mergers at z < 1 are mostly gas-poor, 
while mergers with a higher mass ratio usually involve 
a more gas-rich lower mass galaxy and so minor and in- 
termediate mass ratio mergers are gas-rich until lower 
redshift. This is not necessarily t he case in other models 
(see Figure 7 in |Lotz et ar]|2011l ). 

The merger fraction for the SOB model was found to 
compare well with the observe d fractio n of morphologi- 
cally disturbed galaxies in Jogee et al. (2009) and with 



the major close pair fraction m Lotz et ai 



pUTT] ) (see also 
earlier versions 



|Bell et aL\ 2006a1b| for comparisons with 
of the model). 

5.2. Mass-Metallicity relation 

The normalization and slope of the color-magnitude re- 
lation produced by the model are strongly influenced by 
the M-Z relation, but also depend on the star formation 
histories of the galaxies. A line ar fit to the observed rela- 
tion from Gallazzi et al. ( 2006 ) results in good agreement 
with the slope ot the observed red sequence and a small 
offset in color compared to observations. We adjusted 
the fitted M-Zrelation by 0.004 (see Equation (m)) to 
improve the match at z = 0.9. The mean color of tne red 
sequence is slightly redder and no offset to the M~Z- 
fit is required if we do not incorporate significant bursts 
of star formation but assume only a consta nt st ar for- 
mation history for each galaxy (see Section 5.3). The 
consequences of dry mergers for the colors of massive red 
galaxies depend on the existence of a tilt in the color- 
magnitude relation but are fairly insensitive to the exact 
choice of M-Z relation. 

The scatter in the red sequence is driven mainly by 
age differences. The scatter in all three models is smaller 
than observed. This leaves room to include scatter in 
the M-Z relation. Assuming a log normal relation with 
a wid th of 0.1 dex (Tremonti et al. 2004 Gallazzi et al. 
2006 1 results in slightly larger scatter m ail the model 
red sequences, still well within the range expected from 
observations. The differences between the models are 
evident even when improbably large scatter in the M-Z 
relation is used. 

One could imagine that if all mergers back to the high- 
est redshifts are considered, most dark matter halos are 
likely to have had a major merger, with many of these 
mergers occurring at early times when the masses of the 
galaxies involved were very small. Using a non-evolving 
M-Z relation then implies that the metallicities would 
be very low. With no further star formation, the metal- 
licities of remnant galaxies on the red sequence would 
remain unrealistically low. This turns out not to be a 



concern. We consider mergers taking place after the as- 
sumed formation redshift of z/ = 4. Furthermore, we 
assume that only the final gas-rich major merger on any 
branch of the merger tree moves the remnant onto the 
red sequence. This excludes some very gas rich merg- 
ers occurring at high redshifts, which is a reasonable as- 
sumption as long as the remnant retains sufficient gas to 
remain in the blue cloud or re-accretes gas before a sec- 
ond major merger. It does lead to a concern that similar 
types of mergers may have different effects depending on 
when they occur, or be entirely ignored in one case while 
creating a red sequence galaxy in another case. An alter- 
native approach could be to enforce an upper gas frac- 
tion threshold, above which the remnant is assumed to 



remain in the blue cloud (see Section 5.4 below). These 



conceptual issues largely affect the creation of the model 
red sequence but will have little effect on the differences 
between the models at low redshift. We have shown that 
the red sequence produced by z = 1 under the current set 
of assumptions agrees sufficiently well with observations 
to merit its use as a reference. 

5.3. SF histories 

We have chosen a very simple constant star formation 
history for all galaxie s as the defa ult. This was also the 

f 



approach followed by Kriek et al 
of red sequence color evolution. 



2008| in their models 
here are many other 



options that could be c onsidered. The scatte r we find is 
consistent with that of Ruhland et al. ( 2009) , who used 
an exponentially declining Sl'K in their exploration of 
how the scatter in the OMR is affected by the continu- 
ous addition of quenched blue cloud galaxies. An expo- 
nentially declining SFH may not be a good description 
for all galaxies, however. Lyman break galaxies that are 
forming stars at high levels at redshifts above two may 
be better described by an increasing S FH, for example 



(Lee et al. 2010 Papo vich et al. 2011). Previous work 



on the progenitor bias ( TO.n Dokkum &: l''ranx|2001 ) used 



three parameters to describe the star formation history 
— the start time, end time, and a parameter that de- 
scribes how SF varies within this period. They note that 
the details of how the SFR changes with time are not 
very important; the evolution of an individual galaxy can 
be approximately described as the evolution of an SSP 
that formed at the time giv en by the mean luminosity 



weighted age of the galaxy (van Dokkum et al. 199Sb 



van Dokkum & Franx 2001) 



1 



In our fiducial model, star formation for all galaxies 
begins at a redshift of four. We tested more extreme 
starting times corresponding to 2; = 2 and z = S, as well 
as a variation in which the starting times are randomly 
chosen from a Gaussian distribution of width 0.5 Gyr 
centered on z = 4 and constrained to lie between z = 10 
and z = 1.5. For a more recent starting time {zf — 2) 
the model red sequence is bluer than the observed rela- 
tion by approximately 0.05 mag aX z — 0.9. At z = 0.1, 
the red sequence for Model I is slightly redder than ob- 
served, Model II matches well with the observed relation, 
and Model III remains slightly bluer than observed. The 
luminosity evolution for all three models is faster as all 
the stars have formed more recently. This gives larger 
changes in both color and magnitude between z = 0.9 
and z = 0.1. The relative changes for the three models 
are in the same sense as before but the impact of allowing 
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star formation to continue after z = 1 is lessened, leading 
to a smaller difference between Model II and III. When 
star formation starts at z = 8, the evolution is slightly 
slower for all models but the results are very similar to 
the fiducial case. The results when a range of starting 
times is allowed are almost indistinguishable from those 
with a single starting time. 

If major wet mergers are the catalyst for the quenching 
of star formation, as we have assumed in our models, they 
are likely to be accompanied by a burst of star formation 
that uses up or expels the remaining fuel from the central 
regions of the galaxy. We found that it was necessary to 
convert the significant amounts of cold gas in major wet 
merger progenitor galaxies in the SAM to stars in order 
to produce high enough masses. It is not surprising that 
neglecting the cold gas results in lower normalization, as 
this gas is converted into stars in the full SAM, and the 
MF of the SAM at z = matches well with observa- 
tions. Bursts of star formation have stronger impact at 
higher redshifts, where the galaxies have higher gas frac- 
tions, approximately doubling the mass of galaxies with 
^ 3 X 1O^°M0. We tested other ways of implement- 
ing bursts as well as star formation histories in which 
the same amount of mass is formed with no burst. It has 
been shown that star formation in mergers is only mildly 
enhan ced with respect to normal star-forming galaxies at 



z ^1 { Barton et al. 2007 Robaina et al. 2009 Jogee et al 



2009 1 ; thus we would not expect a large fraction of the 



mass to form in merger-induced bursts of star formation. 
The enhancement in star formation due to mergers in 
the S08 SAM was fou nd to agree well with observations 
(Robaina et al. 2009), suggesting that dramatic bursts 
of star formation are not needed, at least at z < 1. At 
higher redshifts, more intense bursts of star formation 
may occur. In general, the enhancement in SFR caused 
by a merger is dependent on the morphologies, masses 
and gas fractions of the galaxies involved, as well as the 
orbital parameters and the mass ratio. The relative evo- 
lution of the three models and our conclusions remain 
the same whether the gas mass is excluded, included in 
a single burst at the time of the merger or incorporated 
into an extended period of star formation, although the 
colors do depend on this choice. With a constant star 
formation rate and no burst, the colors are redder by 
^ 0.1 mag on average. The color is also dependent on 
the length of the burst; a shorter, more intense burst of 
star formation produces slightly bluer galaxies, but again 
has little impact on our conclusions. Allowing star for- 
mation to continue for some time after the merger may 
be more realistic than the sudden quenching we assume, 
but it would make little difference to the results. 

It is possible to produce a CMR with a similar slope 
and normalization aX z — \ (using the same M-Z rela- 
tion) by assuming that all the stars formed in a single 
burst at fairly high redshift (z ^ 2) rather than over 
an extended period of time. The scatter in the relation 
is caused by differences in the color evolution for dif- 
ferent metallicity populations rather than differences in 
age. We find that the scatter at the massive end of the 
relation in this model is unrealistically small at z = 0.9 
(0.02 mag for A/* > IO^^Mq). Allowing for scatter in the 
Af-Z-relation would compensate for this to some extent. 
Although the scatter in the relation as a whole remains 



approximately constant to low redshifts, the scatter for 
the most massive galaxies actually increases slightly in 
this model (0.04 mag at z = 0.1). The evolution in color 
and luminosity to low redshifts is too rapid compared to 
observations, even if dry mergers are included. An earlier 
formation redshift [zj — A) results in too red a relation 
by z = 1 although the change in color between z = 1 and 
z = is approximately correct. This incompatibility be- 
tween high redshift single-burst passive evolution models 
and the observed c olor evolu tion has been pointed out by 
Kriek et al] ( |2008[ ); [Whitaker et al.| ([2010) and others. 



5.4. The shape of the mass and luminosity Junctions 

The shape of the model LF and MF differs somewhat 
from the observed relations (see Figure [2]). The number 
densities of galaxies at the bright end and faint end in 
the models that include mergers after redshift one (Model 
II and III) ag ree reasonably with the observed number 



densities from Brown et al. ( 2007 ) but there is a shortage 
of intermediate luminosity galaxies. In the models, the 
evolution in the MFs results solely from mergers, but in- 
cludes new additions to the red sequence through major 
mergers, as well as the increase of mass in existing red se- 
quence galaxies. There is very little evolution below the 
knee of the MF, even in the two models where mergers oc- 
cur, suggesting that more migration of blue cloud galax- 
ies onto the red sequence than predicted by this model 
is necessary to match the observed evolution at interme- 
diate masses. Secular processes such as disk instabilities 



bars and clumps (e 


j;., pfstathiou et al. 
CoTe et al..|2000i ,Bo 


1982; 


Mo et al. 


1998 


Combes||2000 


wer et al.|2U0t)p. 



too early. 

The masses of merging galaxies are extracted from the 
SOS SAM, so similar discrepancies in the shape of the 
MF are also expected for the full SAM. The total MF of 
the SAM, as well as the MFs split by color, agree w ell 
with observations from the SDSS (Bell et al. 2003) at 
low redshift. There is very little evolution in the SAM 
MF, however, resulting in increasing differences to the 
observed relations at higher redshifts. This seems to be 
the case in a number of SAMs and is discussed in detail in 



Fontanot et al. (2009). Although the SAM produces a bi- 
modal distribution of galaxies in color-magnitude space 
at z = 0, the distribution also becomes increasingly dif- 
ferent from the observed distribution at higher redshifts. 
Even if the total LF of the SAM matched well with ob- 
servations at z ~ 1, the relative contributions from red 
and blue galaxies would not necessarily agree. The dif- 
ficulties in matching the observed LFs in detail were to 
some extent our motivation for simplifying the model 
and using only the fairly robust merger trees from the 
SAM as the basis of this work. Future work using the 
halo occupation distribution (HOD) to ensure that the 
distribution of galax y masses in the model matches that 
of ob servations (e.g. iHop kins et aL||2008b Conroy et al. 
20071 will allow us to compare quantitative predictions 
with the observed evolution. In this paper we concen- 
trate on the relative changes between different models at 
a fixed mass and space density to demonstrate the effects 
of mergers, rather than making a direct comparison with 
the observed evolution. 
It has been pointed out that bulges are less likely 
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to form in extremely gas-rich mergers (Hopkins et al. 
2009b|c| ). As less massive galaxies have higher gas trac- 
tions, this would inhibit the formation of low-mass red 
sequence galaxies. This effect is redshift dependent be- 
cause galaxies had higher gas fractions in the past. The 
efficiency of bulge formation depends on mass ratio, gas 
fraction and the orbital parameters of the merger. It is 
beyond the scope of this work to include this level of de- 
tail, however we found that applying a simple upper gas 
fraction threshold for spheroid formation strongly sup- 
presses the number of low mass red sequence galaxies 
formed. We find excellent agreement between the ob- 
served and model MFs at the low-mass end at z = 0.1 us- 
ing an upper gas fraction threshold of 80%, above which 
the merger remnants are assumed to stay in the blue 
cloud rather than moving onto the red sequence. The 
MF dX z — 0.9 overshoots the observed distribution at 
low masses, however this is where incompleteness in the 
observations may be an issue. If no upper threshold is 
applied, the LF and MF overshoot the observations at 
low masses. With a threshold of 70%, the mass functions 
agree well at z = 0.9 but more fading of low mass blue 
cloud galaxies without mergers is required after z = 1 
to reproduce the low-redshift mass function. Adjusting 
the threshold has almost no impact on the growth at the 
massive end. 

5.5. Choice of stellar population model and IMF 

We have verified that the implications of mergers on 
the evolution of color and luminosity described above are 
independent of the exact choice of stellar population syn- 
thesis model. Using the models of Maraston ( 2005 ) with 
a Kroupa IMF results in somewhat slower evolution be- 
tween z — 1 and the present day for simple stellar popula- 
tions, and correspondingly slower evolution in both color 
and magnitude in all three models. The color-magnitude 
relation, LF and MF are similar to those shown in Fig- 
ure [2j and the differences in evolution between the mod- 
els at low redshifts are in the same sense a s shown in 
Figure The main difference between the |Bruzual~fc 
CliarlotT |20"03') and 'Maraston' ("2005) models is the treat- 
ment or post-main sequence evolution, particularly TP 
AGB stars. The impact is strongest in the near infrared 
bands and a t high redshifts, as shown by |van der Wei 



et al. (120061 and others, and so we do not expect very 
big differences. 

The choice of IMF affects the evolution of both lu- 
minosity and color. The distribution of stellar masses 
has a stronger impact on the luminosity evolution than 
the color evolution, which is determi ned largely by the 
age of the stellar population, however ( Tinsley[1980 van 
Dokkum||2008). The magnitude of the evolution iii our 



modeis would depend on the IMF, but we expect the 
differences between the models to hold for any of the 
"standard" IMFs. 

5.6. Relation to previous works 



The progenito r bias ([van Dokkum fc Fr anx''1996 



Dokkum et al. 



van 



2000| [van Dokkum I'ranx^ 200f| ^ W 



consequences ot mergers that we have demonstrated here 
have been recognized as an important factor for the evo- 
lution of the red sequence for some time and pointed out 
in numerous papers on the evolution of massive early- 
type galaxies (for e.g.. Brown et al.|2007 Bell et al.|2007 



Kriek et al.]|2008[ [Whitaker et al ||2010[ ). What is new 
in this work is the calculation of tiie magnitude of those 
effects in a hierarchical model, and a comparison of the 
evolution of massive red galaxies with and without merg- 
ers after z ^ 1. 

( |2006a|b l 



Models by Hopkins et al. 
co-evolution of early-type ga. 



considered the 
axies and quasars, assum- 



ing a quasar phase is triggered by merging activity and 
spheroidal galaxies form as remnants of mergers. They 
use the observed quasar luminosity function as a starting 
point and reproduce a number of properties of early-type 
galaxies, including the LF, MF and evolution of the M/L 
ratio. Later models using the HOD approach with the 
simple assumption that major mergers are responsible for 
quenching star formation and forming spheroids also pre- 
dicted a number of properties of early-type galaxies tha t 
match well with observations (Hopkins et al.''2008b"a). 
We have used a simpler model framework but make sim- 
ilar assumptions on the formation of spheroids. Our re- 
sults are in good agreement. 

The combined effects of mass build up through merg- 
ers and stellar evolution that we test here ar e naturally 



included in more complex SAM s (see, e.g 
eraL|[2008l [Bower et al 



2006 



S omervillel 
De Lucia & Blaizot 



irvillej 
2007f 



Monaco et al.||2007p . Many of these models have implic 



itiy or explicitiy shown that massive early-type galaxies 
form hierarchically rather than evolving purely passively 
for many years now. In particular, it has been shown that 
the fairly rare massive early- types in the densest environ- 
ments (brightest cluster galaxies) continue assembling to 
low redshifts, mostly through dry mergers, although the 
stars in these merging components m ay have formed at 
high redshifts ([Aragon- Salamanca et a l. 1998; Dc Lucjal 
fc Blaizot|2007p . In a recent paper De Lucia et al.| ( |201iy 



examine how much mass in bulge galaxies has been con 
tributed by mergers and disk instabilities. Interestingly, 
they find that disk instabilities play the most important 
role at intermediate masses, where the discrepancy be- 
tween the model and observed MFs is greatest in our 
merger-only bulge-formation scenario. 

The LF is often used to probe the mass growth of mas- 
sive early-types by testing whether there is evidence for 
any evoluti on beyond that of an old, passively evolving 



model (e.g. |Wake et al.||2006[ Cimatti et al.||2006[ |Cool 
jeTal. 2008). We have shown that even with mergers the 
LI' can appear to evolve only by the amount expected 
from passive evolution. One of the strongest implica- 
tions of this work is then that the LF alone cannot be 
used to distinguish between the passive and hierarchical 
formation scenarios for massive early-type galaxies. Mass 
measurements that are independent of the assumption of 
passive evolution are required. A recent measurement of 
the evolution of the MF where the LF is converted to 
a MF using the observed evolution of the M/L led to 
a higher estimate of the change in number density of 
massive galaxies that is more consist ent with the p redic 
tions of hierarchical models ( Robain a et al.|[20T0 ) . Mat 



suoka fc Kawara ([2010") also measure a greater change 



in number density using masses determined from spec- 
tral energy distribution (SED) fitting, but other works 
that use similar methods still find very lit tle evolution 
([Ferreras et al.||2009[ [Banerji et"aI1|2010[). Note that 



there is strong model dependence in htting SEDs to ob- 
tain masses, which can lead to large uncertainties (see 
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Gallazzi fc Beli||2009 ). 



6. CONCLUSIONS 

We have used galaxy merger histories from a hierar- 
chical model, in conjunction with simple star formation 
histories and stellar population synthesis models, to ex- 
plore the color and magnitude evolution of red sequence 
galaxies in several scenarios. 

We use the galaxy merger trees from the SOS S AM com- 
bined with stella r evolutionary tracks from the |Bruzual] 
fc Chariot ( 2003 ) model, making simple assumptions on 
the relation ot mass to metallicity and the redshift at 
which galaxy formation begins. We assume that major 
wet mergers are effective at quenching star formation but 
that star formation may also shut down after some pe- 
riod of time, independently of mergers. With this model 
setup, we examine the role of recent mergers and the ef- 
fect of adding younger populations to the red sequence 
at late times. 

A red sequence population with a realistic color- 
magnitude distribution at z = 1 can be formed through 
the passive fading of galaxies that form in a single burst 
relatively recently {zf ~ 2) or through more extended 
star formation, shut down by major mergers at different 
times, in galaxies that began forming earlier (zf ~ 4). 
Galaxies that form all their stars at high redshifts have 
colors that are too red by z = 1 and produce a CMR that 
has smaller scatter than observed, but the slow redden- 
ing they experience thereafter gives approximately the 
observed amplitude of color evolution. 

A population with approximately the correct color at 
z = 1 evolves very rapidly over the last half of cosmic 
history if no further mergers occur and there is no new 
mass added to the red sequence (Model I). This results 
in a color-magnitude relation that is too red at low red- 
shift, with too few bright red galaxies and too much lu- 
minosity evolution at fixed space density compared to 
observations. At a fixed mass, the change in color and 
magnitude predicted by such a model are approximately 
the same as the changes predicted for a purely passively 
evolving population that formed at z/ ~ 2. 

Dry mergers occurring between z = 1 and z — slow 
down the luminosity and color evolution of the early- 
type galaxy population. The mass (and light) added 
through mergers counteracts the fading of the luminosi- 
ties expected from stellar evolution when no new stars 
are forming. The result is a slightly bluer CMR at z = 
that is closer to the observed relation, and more realistic 
evolution of the LF and MF (Model II). By allowing star 
formation to continue in some galaxies after z — 1, the 
changes in a population's color and luminosity are re- 
duced even further (Model III). This late morphological 



transforination is known as progenitor bias ( van Dokkum 
fc Franx||2001[ ). Both the CMR and the LF evolution ot 
this model are consistent with the observed evolution. 
The smaller changes in color and magnitude at a fixed 
mass replicate the behavior of a passively evolving pop- 
ulation that formed in a single burst at high redshift 
(z/=3-5). 

At the bright end of the red sequence, the seemingly 
slow passive evolution of luminosity and color is equally 
well interpreted as the passive evolution of an ancient 
population, or as the cosmologically motivated hierar- 
chical growth of an evolving population. Although they 



predict very similar evolution, the implications of these 
scenarios for the growth of stellar mass are very differ- 
ent. Diagnostics other than the LF, such as the num- 
ber density evolution estimated using dynamical M/L 
information, the growth rate determined from close-pair 
statistics and the evolution of the clustering strength of 
early-type galaxies lend support to the hierarchical pic- 
ture. 
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